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Abstract 



00 ■ We develop a systematic analytic approach to aging effects in quantum 

disordered systems in contact with an environment. Within the closed-time 



path-integral formalism we include dissipation by coupling the system to a 



set of independent harmonic oscillators that mimic a quantum thermal bath. 



After integrating over the bath variables and averaging over disorder we ob- 
tain an effective action that determines the real-time dynamics of the system. 
The classical limit yields the Martin-Siggia-Rose generating functional asso- 
ciated to a colored noise. We apply this general formalism to a prototype 
model related to the p spin-glass. We show that the model has a dynamic 
phase transition separating the paramagnetic from the spin-glass phase and 
that quantum fluctuations depress the transition temperature until a quan- 
tum critical point is reached. We show that the dynamics in the paramagnetic 
phase is stationary but presents an interesting crossover from a region con- 
5h ' trolled by the classical critical point to another one controlled by the quantum 

critical point. The most characteristic property of the dynamics in a glassy 
phase, namely aging, survives the quantum fluctuations. In the sub-critical 
region the quantum fluctuation-dissipation theorem is modified in a way that 
is consistent with the notion of effective temperatures introduced for the clas- 
sical case. We discuss these results in connection with recent experiments in 
dipolar quantum spin-glasses and the relevance of the effective temperatures 
with respect to the understanding of the low temperature dynamics. 
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I. INTRODUCTION 



In the last few years, great progress has been made towards the understanding of the 
out of equilibrium dynamics of classical systems with slow dynamics and aging effects. Ex- 
periments in several glassy materials such as polymer glasses,! spin-glasses,i orientational 
glasses,!! simple liquids like glycerolB and gels,! have shown that all these systems share 
an extremely slow dynamics at low temperatures. Although some aspects of the dynamic 
evolution as the precise scaling laws or the dependence on cooling rates vary from system 
to system, all the examples mentioned above are characterized by the existence of a non- 
equilibrium low temperature phase with aging effects. This means that measurements show 
a strong dependence upon the time elapsed since its preparation by a thermal quench or an- 
nealing from a high temperature phase. Needless to say, the standard equilibrium approach 
is not suited to describe such a phenomenology. 

In order to describe the spin-glass number of theoretical methods have been 

exploited, ranging from scaling arguments based on coarsening ideas;! phenomenological 
models founded on particular assumptions on the structure of phase space; analytical 
solutions to models with infinite-range interactions: and numerical simulations of finite 
and infinite dimensional models.^ All of these approaches have succeeded in capturing some 
of the effects observed experimentally. 

Importantly enough, the solution of simplified models with infinite range interactions 
has been fruitful for several reasons: 

(i) It has shown that indeed a large family of such models do capture aging effects that 
resemble, at least qualitatively, the experimental observations. 

(ii) It has provided a general framework to analyze the non-equilibrium dynamics with 
new predictions, like scaling laws for the time-delayed correlation function and asymptotic 
violations of the fluctuation-dissipation theorem, that are not restricted to the infinite di- 
mensional case. Some of these have been confirmed by numerical studies of more realistic 
modelsEl and will be probed experimentally soon.0 

(iii) It has allowed to extend the analytical connection^ between mean-field theories of 
disordered systems and the mode-coupling theory of super-cooled liquids^ to the low tem- 
perature glassy phase. This is most welcome since many of the similarities observed in 
the behavior of real glasses of very different origin can now be understood theoretically. 

(iv) Problems of a more general nature such as the dynamics of manifolds in a random 
quenched environment, relevant to the description of e.g. the pinning of the flux lattice by 
disorder in dirty high T c superconductors,!! can be studied in a similar way.0 The depen- 
dence on the internal dimension of the manifold can be included in the formalism and the 
analysis yields, among other results, dynamic scaling laws for the displacements with aging 
effects. 



The studies we have mentioned so far are mostly concerned with the analysis of classical 
systems. This treatment is in general justified since, in many cases, the critical temperature 
at which the transition to the "glassy phase" occurs is sufficiently high as to make quantum 
effects irrelevant. Nevertheless, in many cases of practical interest, the critical temperature 
can be lowered by tuning an external parameter. Experimental realizations of this situation 
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are given by insulating magnetic materials^ such as LiHo x Yi_ x F4 and randomly mixed 
hydrogen-bonded ferro-antiferroelectric crystals of the Rbi_ x (NH^H^PC^ type (known as 
RADP)B 

Theoretical studies addressing the effects of quantum fluctuations in disordered media 
have been considered over the last twenty years. A variety of techniaues including replica 
theoryEM renormalization groupie and Monte Carlo simulationscatill have been employed 
in order to understand the static properties and critical behavior of quantum systems such 
as quantum rotors!!, transverse field Ising models in finite dimensional! and mean-field 
systems like the Sherrington-Kirkpatrick model in a transverse field.il~I! 

From the study of classical glassy systems described in the first paragraph, one knows 
that at low temperatures their dynamics becomes so slow that the systems are not able to 
reach thermal equilibrium at any experimentally accessible time. The same occurs regarding 
numerical simulations of large systems. It is then plausible that even at very low tempera- 
tures, say T = 0, with only quantum fluctuations driving the dynamics, the materials and/or 
models will need an unaffordable long time to reach equilibrium and that the relevant evo- 
lution will be characterized by non-equilibrium effects. Indeed, it is well-known that certain 
glasses in the limit of very low temperatures have aging effectsBll and hence the approaches 
based on the assumption of equilibration have a restricted domain of application. Instead, 
one should necessarily start from a quantum dynamic description in order to obtain sensible 
information about the systems. To the best of our knowledge this problem has not been 
theoretically tackled, from a microscopic point of view, yet. 

The question then arises as to how do quantum fluctuations modify the real time dy- 
namics of glassy systems. In this paper, we treat the case of a disordered model in contact 
with a thermal quantum bath. 

Our first aim is to present a general framework that allows us to study the real time 
dynamics of a system with quantum and thermal fluctuations starting from an arbitrary 
initial condition. In the case of classical models this is most easily done within the Langevin 
approach. One generally starts by postulating that the dynamics of the microscopic variables 
is given by a stochastic Langevin equation. The noise and friction terms account for the 
bath-system interaction. A path-integral formulation of the generating functional, known as 
the Martin-Siggia-Rose formalism^ (MSR), yields an elegant and simple formulation which 
is particularly well suited for problems with disorder. However, the generalization of this 
procedure to the quantum case is far from being straightforward. The Langevin equation 
does not seem to be a good starting point since, being phenomenological, it cannot be 
quantized in an obvious way.0S In other words, it is not clear how to make the Langevin 
approach compatible with the intrinsic microscopic dynamics of a quantum system. 

The approach we use here results from the combination of different techniques which 
have been intensively discussed in the literature and successfully applied to a wide range 
of phenomena. Essentially, it consists of the Schwinger and KeldyshU closed-time path- 
integral formalism applied to a system in contact with a thermal bath, combined with the 
integration of the bath variables leading to the Feynman and Vernon influence functional. §A 
subsequent integration over the disorder gives the quantum analogue of the MSR generating 
functional for disordered systems. 

Though the formalism allows us to consider general cooling procedures or/and various 
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temperature variations, we shall not study involved thermal histories here but concentrate 
instead on a quench from a configuration associated to the disordered phase. 



For the sake of definiteness we focus on a system whose classical counterpart, the so- 
called "p-spin spherical spin-glass model", has been studied in great detail. This model is 
described by a potential energy 

N 

V[tr,J] = - E J i i ...iv a h--- (T if ( L1 ) 

il<...<i p 

with Ji lt ...,i p quenched independent Gaussian variables and p > 2 a parameter. 

Let us attempt a short review of the behavior of the classical p-spin models. 

In its Ising version it was first introduced by Derrida in relation with the random en- 
ergy model.0 Several studies of its static properties, in particular using replicas, followed.0 
Khkpatrick and ThirumalaB studied the Equilibrium dynamics, a la Sompolinskypi@ and 
remarked that the dynamic equations, in the high temperature phase, are those appearing 
as F p _i models in the mode-coupling theory of super-cooled liquids developed to describe 
the glass transition.^! 

The spherical version, where a time- dependent Lagrange multiplier is introduced to 
enforce the spherical constraint J2iLi °f = N, was first studied by Jones, Kosterlitz and 
ThoulessB for p = 2. Its statics is fully solvable either by using replicas - a replica symmet- 
ric ansatz is exact - or through a direct calculation. Its dynamics^i@ is also fully solvable 
and it renders explicit the fact that the model is formally relatedffl to the 0(N) model 
in the large N limit and dimension d = 3. In the classical limit, it is then a toy model 
for ferromagnetic domain-growth. In the quantum case one feels tempted to relate it to 
coarsening in quantum systems. 

The analysis of the spherical model for p > 3 is due to Crisanti and Sommers,0 who 
also studied the equilibrium dynamics a la Sompolinsky in collaboration with Horner.@ In 
the classical limit and at the static level, this model is the prototype spin-glass solved by a 
Parisi one-step replica symmetric ansatz. Still at the classical level and in the large N limit, 
one derives exact evolution equations that exhibit a first order - discontinuous - dynamic 
transition at TV Above T^, as in the Ising case, one recovers the mode-coupling equation 
for model Below Td it is well-known that the dynamic approach used in Ref. |@] is 

incorrect .0^0 The dynamics in the low-temperature phase was solved in Ref. JIT] showing 



that the dynamics of mean-field spin-glasses strikingly resembles that of real systems. In 
the spin-glass phase the evolution is intrinsically out of equilibrium and the systems age in 
similar way to the one observed in experiments. 

In order to define the model at the quantum level, we need to specify the dynamics of 
the variables a . For simplicity, we consider a kinetic term of the form 

K = yf> 2 (1.2) 
z i=i 

where the dot denotes a time-derivative. Calling z the Lagrange multiplier enforcing the 
spherical constraint, the model we here consider is then described by the action 
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S S [CT,J] = Jdt 



m 



■2 



N 



N 

i\<...<i p 



;i.3) 



Our main goals are to detect a dynamic transition from a disordered phase to an ordered 
glassy phase and to characterize the dynamics in these phases. With this purpose, we pay 
particular attention to the departure of correlation and response functions from station- 
arity and, we search for one of the hallmarks of a glassy phase, namely, violations of the 
fluctuation-dissipation theorem. We also intend to discuss our results in view of previous 
and possible experimental tests. 

The paper is organized as follows. In Section || we discuss the closed-time path-integral 
formalism of Keldysh and Schwinger and the coupling of the quantum system to a quantum 
heat bath. We obtain a generating functional with an effective action that takes into account 
the effect of the coupling to the bath. In three Appendices we give details of the Keldysh 
Schwinger formalism (Appendix |S|), we recall the quantum FDT (Appendix |B|) and we 
discuss the classical limit (Appendix [C]). Section |TTT| is devoted to the derivation of the 
dynamic equations through a saddle-point approximation of the generating functional. The 
method is similar to the one used in the classical case when manipulating the Martin- Siggia- 
Rose generating functional. The general scenario for the real-time dynamics of this system 
is discussed in Section [IV]. It is inspired by the behavior of the classical counterpart but 
some generalizations are necessary. In Section [V] the equations and their solution in the 
paramagnetic phase are studied with special emphasis set on the analysis of the out of 
phase susceptibility close to the transition line. In Section [VI| the dynamic equations in 
the glassy phase are studied. The non-equilibrium dynamics is shown to survive quantum 
fluctuations in a finite region of the phase diagram. Further details of the approach to solve 
the equations are given in Appendix [B]. The modification of FDT needed to construct an 
ansatz to solve these equations and its relation to effective temperatures are discussed in 
Section |V11| . Finally, in Section |V11I| we present a summary of our results, our conclusions 
and some ideas as to how to continue investigating along the lines here proposed. A short 
account of some of these results has appeared in Ref. Wl . 



II. CLOSED-TIME PATH-INTEGRAL FORMALISM 



The closed-time path-integral (CTP) formalism of Keldysh and Schwinger^l for treating 
the non-equilibrium dynamics of a quantum system has been reviewed several times in the 
literature.E3tia In order to set the notation we briefly present it in Section |11 A| . We then 
consider the situation in which the system under study is coupled to a quantum thermal 
bath and derive the effective action following Feynman and VernonS (Section [II B|) . In 
Section |II C| we finally address the case of a disordered system and we obtain the effective 
action from which the dynamic equations are derived in Section ID. 
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A. General framework 



The basic ingredient of this formulation is the "in-in" generating functional from which, 
after derivation with respect to the external sources 4 + , the different Green functions 
can be obtained. Let us consider a model with degrees of freedom described by a field 
<p(t) = (0i(t), . . . , 0jv(O) m the Heisenberg picture. Denoting p(t Q ) the density matrix at 
the initial time t (which from now on we set to zero), the generating functional is defined 
as 

Z\e,C] = Tr [T* exp (-1 jT dtC(t)<t>(t)) Te W Q dt )*(*)) p(0)] (2.1) 

where £ + {t) and are the N- vector external sources and a scalar product between N 

vectors is assumed. The symbols T and T* are the time and anti-time ordering operators: 

(2.2) 
(2.3) 



TiAitjBfo)) = 9(h - t 2 )A(t 1 )B(t 2 ) + 9{t 2 - tjB^Afa) 
T*(A(h)B(t 2 )) = 6(t 2 - tJAitjBih) + e(h - t 2 )B(t 2 )A(t 1 ) 

Defining ensemble averages in the usual way 

Tr(A(t)p(0)) 



(A(t)) = 



one obtains 



S 2 lnZ 



4=o 



Trp(0) 



which is directly related to the Feynman causal propagator. In the same way, 



1 

2 

2 



c ij (t,t') = -(<t> i (t)<j> j (t') + (*)> 

+ 



\nZ 



(2.4) 
(2.5) 

(2.6) 
(2.7) 



where Cij(t,t') is the symmetrized correlation function (Cij(t,f) = Cji(t',t)). 

The response function Rij(t, t') is defined as the variation of the averaged field {(pi) with 
respect to the "strength" fj of a perturbation that modifies the potential energy according 
to V -> V - f<f>: 



f=o 



In linear-response theory, it can be expressed in terms of the averaged commutator: 



R iJ (t,t') = ^-6(t-t')([ ( j )t (t),Mt > )]) 



% 
h 



\nZ 



(2.8) 

(2.9) 
(2.10) 
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One of the main advantages of this formalism is that the generating functional (|2.1|) 
admits a path-integral representation which is a slight modification of the usual "in-out" 
Feynman path integral. This will allow us to make formal manipulations similar to the ones 
performed in the study of the dynamics of classical systems via the Martin-Siggia-Rose@ 
(MSR) formalism. In a short-hand notation, the generating functional reads 

Z[t,C] = J Vcf> + Vct>- exp [~ - S[4T] + / dtt(t)cf> + (t) - J dtC(t)<t>-(t)) 

x(0+|p(O)|<r> (2.11) 

where S[<f)] is the action of the system and (0 + |p(O)|</>~) stands for the matrix element of 
the density matrix which has support only at t = (see Appendix [A] for a proof of this 
equation). Hereafter, we omit the limits of the time-integrals that, unless otherwise stated, 
go from t D = to oo. 

In the following it will be useful to write the correlation and response in terms of the 
fields <p + , <f> : 

^i(M') = i (<Pt(t)cf>j(t') + ^{t')cj>T(t)) , (2.12) 

%(M') = \ (tfit) (#(0 - 07(0)) • (2-13) 

Notice that as a result of having two external sources (one related to the evolution 
with time ordering and the other with anti-time ordering) a doubling of degree of freedom 
is necessary when writing the path-integral. As we explicitly show in Appendix 0, the 
doubling of degrees of freedom characteristic of this description is intimately linked to the 
introduction of Lagrange multipliers in the MSR description of the classical model. 

The time-integration in Eq. (|2.11| ) can be interpreted as being closed, going forward from 



to = to t = oo and then backwards from t = oo to t = 0. This motivates the name of the 
method. 



B. Coupling to a heat bath 

Up to now, the discussion has been completely general and it applies to systems with 
an arbitrary number of degrees of freedom. We now consider the coupling of the system of 
interest to a thermal quantum bath, assumed to be in equilibrium. A convenient way to 
deal with the interaction between the system and the bath in the path-integral formalism 
is due to Feynman and Vernon.0 In this approach one starts by describing the system 
plus bath by = (er, v a ) with cr = (ax, ... , ex at) denoting the variables of the system and 
v a = (v", . . . , v%), a — 1, ... , Nf,, denoting the variables of the bath. The total action is 

S[(f>] = S s [<t] + S b [v a ] + S sb [a, v a ] (2.14) 

where S s is the action characterizing the system (and eventually depending upon disorder), 
Sf, is the action for the bath, and S s b contains the system-bath interaction terms. For 
definiteness we assume that our system interacts linearly with the bath which we model as 
a set of independent harmonic oscillators, 
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1 Nb / X ^ 

S[4>] =S s [ct] + -J: M a (yl - col - E Cav a a- . (2.15) 



a=l a=l 



M a are their masses, u a their frequencies and C a the coupling constants between the os- 
cillators and the system. Of course, more general couplings between bath and system are 
possible but this is the simplest choice that allows us to go a bit farther analytically. If we 
further assume that the initial density matrix p(0) is factorizable as 

p(0) = p 8 (0) x p b (0) , (2.16) 

where p&(0) is the Boltzman distribution for the bath variables at equilibrium at a tempera- 
ture ksT = 1/(3, one can integrate out the bath variables, obtaining in this way an effective 
thermal action St[(t + , er~], that enters the Feynman- Vernon influence functional: 



S t [(t + , <r-]=-[dtJ dt' (tr+(t) - *-(t)) V {t - t') (*+(t') + a- (If) 

+*JdtJ dt' (<r+(t) - a-(t)) v(t - t') (a + (t') - tr-(O) • (2.17) 

The noise and dissipative kernels v and r\ are given by@HS 

v(t - t') = J™ dujI(uo) coth (^P?uv \ cos(cj(t - t')) , (2.18) 

poo 

rj(t - t') = -6(t - t') du sin(w(* - if)) . (2.19) 



In these equations, I (to) is the spectral density of the bath: 

I(u,) = ^S( U - Ua )-f-. (2.20) 

a= l ZM a U a 

The system-bath interaction manifests via the appearance of two quadratic and non-local 
terms ( |2.17|) in the effective action. As it can be most easily seen in the classical limit 
presented in Appendix |C] they can be related to the noise and dissipative terms of the 
associated Langevin formulation of the classical dynamics. 

In addition, it is easy to see that the kernels v and rj are related to the correlation 
and response functions of the bath. In fact 4r](t — t') is the response of the bath while 
—2hu(t — if) is the correlation of the bath. They are related in a way dictated by the 
fluctuation-dissipation theorem (FDT): 

Mu) = I hm / — tanh ( 2hu(u') . (2.21) 

(See Appendix [B| for the derivation and properties of the quantum FDT.) This equation 
is the quantum counterpart of the fluctuation-dissipation relation of classical relaxational 
dynamics, that fixes the dissipative coefficient in terms of the strength of the noise correlation 
of the Langevin equation. It expresses the fact that the bath is assumed to be in equilibrium 
at the initial time t Q = and that it stays in equilibrium at all subsequent times t > 0. 
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It is important to notice that the assumption of thermal equilibrium of the bath does 
not imply that the system - that is in contact with it - is in equilibrium at the initial time 
nor that it will reach equilibrium at a given equilibration time t EQ . It is known that "mean- 
field" classical disordered models with relaxational dynamics - the p spin-glass model is an 
example — do not reach equilibrium at any time if the limit of N going to infinity is taken 
at the outset of the calculation.0 It is one of the main aims of this paper to show that 
this result carries through to the quantum case. As a consequence of the out of equilibrium 
behavior, we shall find that the FDT does not hold between response and correlation of the 
system in a sub-critical region of the phase diagram. 

The integrals in the rj and v definitions are, in general, ill-defined (divergent) and might 
need to be regularized. A way to regularize themE3 is to observe that one expects I(co) to 
decrease to zero for large ui and introduce an explicit cut-off: 



I{u) ~ for to > A . (2.22) 

By modifying I (to) in this way, we regularize the kernels 77 and v at the same time preserving 
the FDT relation for the bath variables; in other words, we do not break the assumption of 
equilibrium of the bath. 

Different environments are characterized by different behaviors of I (to) such that for 
uj < A: 

a = 1 Ohmic , 
I (to) ~ uj a and { a < 1 Subohmic , 



a > 1 Superohmic 



A typical example isI^S 



m = ^ *> (jT' 1 exp ( ] (2.23) 



7T 




where M^ is a constant that plays the role of a friction coefficient. 

Throughout this paper, we concentrate in the Ohmic case. The kernel 77 is 

"< ( - ( '> = ^ -'^4^) ( i + ( A(«-f))- ) ' (2 ' 24) 

In the limit where the cut-off tends to infinity it becomes 

lim 7](t - t') = 6(t - t') Mj 5'(t - t') . (2.25) 

A^oo 

Of particular importance is the zero temperature limit of the kernel v: 

T^O V ^ 7T (1 + (A(t-t')) 2 ) 2 

One sees that the quantum fluctuations yield a non-trivial kernel even in the absence of 
thermal fluctuations. In the zero temperature limit the characteristic time goes to infinity 
and the long-time decay of v becomes power-law: 



9 



lim limz/(t-t') = -^(t-t')- 2 • (2.27) 
In the classical limit, for any temperature, v becomes 

lim 2ni/(t - t') = — . 2.28 

R-»0 V ^ 7T l+A 2 (t-t') 2 

Note that for a finite cut-off A, the classical limit of the kernel v also has a rather slow, 
power-law decay at large time differences. If, next, the cut-off is sent to infinity one recovers 
a delta function 



lim lim 2hv(t - t') = AM^ k B T 5(t - t') (2.29) 

A^oo Ti— >0 

characteristic of e.g. a white thermal noise. The classical limit of these terms is further 
discussed in Appendix where we show that the dynamic equations of the corresponding 
classical model correspond to a friction coefficient T^ 1 = 2M / y . 

The coupling to the quantum thermal bath introduces three time-scales into the problem. 
First, the inverse cut-off 1/A characterizes the memory of the bath. Second, the coupling 
strength M^ Q yields the inverse relaxational characteristic time. Third, (f3Ti)~ l accounts for 
the relative importance of quantum to thermal effects. 



C. CTP for disordered systems 

The CTP formalism results to be well suited for problems with disorder.il Let us suppose 
that the system under consideration is described by the action, 

S S [(T, J] = S [(t} - J dtV[tr, J] (2.30) 

where S is the disordered independent part of the action and V[a, J] is the potential energy 
that depends upon random couplings collectively denoted by J. As usual in disordered 
systems, we concentrate on quantities that are averaged over the disorder distribution P[J}. 

We wish to consider the real-time dynamics of the system starting from a random initial 
condition at time t = when it is set in contact with the environment (that is itself in 
equilibrium at a constant temperature T). Given that the initial condition is chosen to be 
random, it is not correlated with the disorder. Therefore, the density operator p(0) does not 
depend upon disorder and the generating functional without sources 

Z\t = 0,r =0,J] =Tr[p(0)] (2.31) 

is also independent of disorder. 

This property is equivalent to the independence of the classical Martin-Siggia-Rose 
(MSR) generating functional^ without sources upon disorder.!^ As in the classical case, 
it allows us to write dynamic equations for random initial conditions without having to 
compute the average over disorder of lnZ[£ + ,£~, J] and hence without resorting to the use 
of replicas. We are then interested in the averaged generating functional: 
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Z[£ + ^-}= J dJP[J) -zfc+r, J), 

from which any averaged operator can be computed as 

1 6Z[£ + ,£-,J] 



(2.32) 



£[0,0, J] 



(2.33) 



cr. Here and in what follows, the overline 



with all sources set to zero and er + = cr 
represents an average over the disorder. 

In many cases of interest, the integration over the disorder can be performed explicitly, 



J dJP[J] exp ~ J dt (V[tr + , J] - V[cr-, J] 



exp 



h 



V D [cr + ,cr- 



(2.34) 



and it introduces a non-linear interaction term that is, in general, non-local both in time 
and in the "space" indices 

A rather general way of introducing randomness is to consider a model with a Gaussian 
random potential energy term V[cr, J] depending upon the variables cr = (01, . . . , <jn) and 
correlated as 



V[cr,J)V[cr', J] = --V 



N ((cr 



N 



(2.35) 



The variables cr may represent spins leading to a spin-glass model or a position in an iV- 
dimensional space becoming the problem of the dynamics and/or diffusion of a particle in 
a random potential. Since we are interested in following the real time dynamics of such 
systems, cr will be a time-dependent iV-vector. 

For the sake of defmiteness, we shall solve a particular realization of the random potential 
that in the classical case defines the so-called p spin-glass:0 



N 



V[cr,J] 



(2.36) 



ii<...<i„ 



p is a parameter, p > 2. The interaction strengths are quenched independent random 
variables with a Gaussian distribution 



P[J] 



\ J 2 irp\ 



exp 



N p-i 



E (Ju 



(4 



J 2 p\ 

2Np- 



(2.37) 



This corresponds to a random potential correlation V(x) = —1/2(1 — x/2) p and it leads to 
an effective, non-local in time, interaction 



V D [cr + ,cr- 



J 2 Ni 
Ah 



dt / dt' 



1 \ p 

j^cr-{t)cr-{t') 



-[^cr+{t)cr-{t') l: 



1 W 
-^cr-(t)cr + (t>) 



(2.38) 
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In its spherical version,0 <7j, i — 1, . . . , N are continuous dynamic variables —\TN < Oi < 
N, Vi, that satisfy the global spherical constraint 



N 



N 



<r\t) 



N 



(2.39) 



i=l 



at each instant. We enforce this constraint by introducing a time-dependent Lagrange 
multiplier z(t). 

The disordered-independent part of the action is 



S [cr, J] = J dt 



•2 



N 



(2.40) 



As explained in the Introduction, for simplicity we have chosen a kinetic term that dictates 
a microscopic dynamics with second derivatives in time. Other choices, for example leading 
to first time-derivatives, are also possible. The derivation of the dynamic equations that we 
explain in Section |TXT| applies with only minor modifications. 



III. MEAN-FIELD EQUATIONS 



We now have all the elements to derive the equations which describe the dynamic evolu- 
tion of a disordered quantum system. According to the discussion in Section ||, the system 
is described by the generating functional 



s K 



J dt$, + (t)(T + (t)- J dtC{t)<T-(t) 



with 



S EFF [cr + ,(T ] = S [cr + ] - S [a- ] + S T [cr + ,(T ]-V D [cr + ,a' 



(3.1) 



(3.2) 



We remind that 5*o is the disorder independent part of the action given by Eq. ( j2.40| ), S T 
accounts for the the system-bath interaction, Eq. ( |2.17 ), and V u is the effective potential 
which arises as a result of the integration over the disorder and is given by Eq. ( [2.381 ) • 

The quantum dynamic equations of motion follow from similar steps to those usually 
used0 to obtain the classical equations of motion in the Martin-Siggia-Rose formalism. We 
introduce macroscopic order parameters and derive, through a saddle-point point approxi- 
mation of the KS generating functional (that becomes exact when N — > oo), the dynamic 
equations of motion. Since we are interested in considering the dynamics of the system for 
large but finite times with respect to N, out of equilibrium effects are expected. We shall 
therefore make no assumption about the time- dependence of the order-parameters. 
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A. Dynamic order parameters 



The quadratic terms in the action can be condensed into one term by introducing the 
operator 



Op(t,t) 



P ++(M') P +-(M') 
P - + (M') Op" (t,f) 



{Op^(M')} 



O v ++(t,t') = (md 2 t + z + (t)) 5(t - t') - 2iu(t - t') 
Op + ~(i, t') = 2r)(t - t') + 2iv(t - f) 
Op~ + (t, t') = -2r](t - t') + 2iv(t - t') 
O v —(t,t') = -(md 2 +z-(t))6(t-t') -2iv(t-t') 

in such a way that 



5' F , PP fcr + , (T 



~\j dt J dt'(T a (t) Ov aP {tit')<r p {t') - V D [(T + } a- 

where Greek indices label a = +, — and the sum convention is assumed. 
Introducing the identity 

1 = jf[DQ^ 5 (l ff a (0c^(O - Q Q/? (M')) 
oc J ]jDQ a/3 D\ a ?exp ( -± A Q/3 (<T a (t)^(t') - NQ aP (t,t' 



into the generating functional, the full action can be rewritten as 



S EFF [er + ,cr 1 



\J dt J dt 1 a a (t) (Op a/3 (M') + X a(3 (t,t')) cr^t') 



+ 



N 
~2 

iJ 2 N 



dt j dt'\ aP (t,t') Q af3 {t,t') + — I dt (z + (t) - z-(t)) 



Ah 



dt / dt' 



(Q ++ (t,t')) P +(Q-(t,t')) P 
Q + -(t,t')) P -(Q-+(t,t')y 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



The stationary-point values of the order parameters Q a ^(t,t') with a, [3 — +, — are related 
to the "physical" correlations and responses defined in Eqs. ( |2.12p and fl2.13| ) as follows 

ih 
~2 

ih 

~2 

ih 

~2 

ih 



{(T+(t)(T- 


-(f)) 


(<T+(t)<T- 


-(f)) 




-(f)) 



NQ—{t,t') = (cr~(t)<T~(t')) = N \C(t,t') + — (R(t, t') + R(t', £)) 



(3.7) 
(3.8) 
(3.9) 
(3.10) 
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with 



1 N 

NC(t,t') = J2(^t(t)cT-(t')+CT-(t)at(t')) 
A i=i 

i N 



NR(t,t>) = £ (<rtV)-°7W)) ■ (3-11) 
11 i=i 

It is easy to check that these functions satisfy the identity [] 

Q ++ (t, t') + Q—(t, f) - Q + -(t, £') - Q~+(t, = 0. (3.12) 

The functional integration over <rf(t) and <J^(t) is now quadratic and can be performed. 
This amounts to replace the quadratic term in i/hS EFF by 

~llJ dt J ^'^gQoP^M') + ~A a/3 (M')) ■ ( 3 - 13 ) 



B. Saddle-point evaluation 

At this stage, all terms in the action depend upon the "macroscopic" quantities \ a @, Q a/3 
and z a and are proportional to N. Since it is easier to write the equations in matrix notation, 
we encode A Q/3 and Q a @ in two matrices 

and we define 

F\o](tt')-( (^(MOr 1 -(Q + -(t,?)y-i\ 

We denote with a cross the standard operational product 

a a ^ _ ( JdtrA+r(t,r)B*W,ir) Idt"A+^(t,t")B^(t",t') \ 
™ 1 ' J \jdt"A-"i(t,t")B"(+{t",t') fdt" A-*i(t,ti')Bn-(if',l>) ) 

where a sum over 7 is assumed. The saddle-point with respect to \ al3 (t,t') yields 

C(t, t') = - Q-\t, t') - Op(t, . (3. 15) 

1 



x At the classical level and in the Martin-Siggia-Rose language, this identity reduces to the con- 
dition (is(t)is(t')) = for all pairs of times t,t'. In the classical case this is a requirement to be 
imposed upon the saddle-point two-point functions in order to preserve causality. In the quantum 
case, this condition emerges as a property of the Green functions. 
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The matrix and time-operator inverse of Q is denoted Q 1 . The saddle-point equation with 
respect to Q a ^(t,t') yields 

C(t,t') = --^F[Q](t,t'). (3.16) 
Equations Q3.15| ) and ( |3.1(j| ) imply, in a compact matrix and time-operator notation, 



t2 

LOp®Q = I-P F[Q] ® Q , (3.17) 
where I is the identity: I af3 (t,t') = 5 a,3 5(t — if). The saddle-point with respect to z a yields 

l = (0p + £)^(t,t) = iQ++(t,t), (3.18) 
% - = (Op + C)zl(t,t) = ±Q—{t,t) (3.19) 

and these equations lead, as expected, to the spherical constraint. 

In the limit N — * oo one could also proceed as in Refs. |I2|,|6"^| and write the full action 
S'eff in terms of a single variable. This is at the expense of modifying the thermal kernel 
and the interaction term in a self-consistent way, through the introduction of terms arising 
from the non-linear interactions (the vertex D and the self-energy E, respectively). 

This procedure is not of particular usefulness for the analysis of the model we treat in 
this paper since the single variable effective action is Gaussian. It does however become 
useful for dealing with quantum models whose single-spin effective action has higher order 
interaction terms. An example is the quantum Sherrington-Kirkpatrick model. One could 
then envisage to derive an effective quantum Langevin equation for the single variablei! and 
study this eauation with an adequate numerical algorithm as the one developed by Eissfeller 
and Opper.a 



C. Dynamical equations for correlation and response 



The dynamic equations for the auto-correlation and response follow from the set of equa- 
tions ( |3.14 )-( ]3.17| ) and the definitions of the dynamic order parameters given in Eqs. ( |3.7| )- 
( ft.lOQ . More precisely, the equation of motion for the response function follows from the 
subtraction of the ++ and H — components of Eq. ( |3.17| ): 



(md 2 t + z + (t)) R(t, if) + 4 /* dt" r]{t - t") R(t", t') 



8(t - 1') 



2ih jo 



dt" [(Q ++ (M")) P_1 - (Q+^MT" 1 )! R(t",t') 



(3.20) 



and the equation of motion for the auto-correlation function follows from the addition of 
the H — and — h components of Eq. (|3.17|) : 
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md 2 t + ~ (z + (t) + «-(*))) C(M') + | - *"(*)) - 

-2% / dt"z/(t - t")R(t', t") + 4 / dt" r)(t - t") C(t", t') 
Jo Jo 



P J 2 



2h 



dt" Im [(Q ++ (t,t")y- 1 Q + -(t",t') - {Q + -{t,t")) p - 1 Q-'{t",t')\ . (3.21) 

Written in this way, Eq. ( |3.21| ) is complex. Its imaginary part yields 

z(t) = z + (t) = z~(t) . (3.22) 

Moreover, since the response is causal, products of advanced R(t, t") and retarded R(t", t') 
responses vanish identically for all t, t": 

R(t,t")R(t",t) = Vt,t" (3.23) 

and one can show that for any integer k > and any constants Ci, c 2 

[C(t, t') + Cl R(t, t') + c 2 R(t', t)] k = [C{t, t') + Cl R(t, t')] k + [C(t, t') + c 2 R(t', t)] k - C k (t, t') . 

(3.24) 

Using this property one has 



(Q ++ (t,0) P-1 " (Q + -(t,t")) p - 1 = Mm 



C{t,t")-^R{t,t") 



p-i 



(3.25) 



and 



Im 



2C(t",t') Im 



-fti2(t',t")Re 



p-i 



We can identify the self-energy £ and the vertex D as 



E(M') +4r/(t-t') = £(t,t') 



p-i 



pJ 2 



t') - 2hu(t - t') = D(t, t') = 



Im C(t,t') - —R{t,t') 

ih 

C(t,t')--(R(t,t')+R(t',t)) 



p-i 



(3.26) 

(3.27) 
(3.28) 



For t > t', they can be encoded in a single complex equation: 

D(t, t>) + | £(t, = ^ (c(t, + | fl(t, o 



p-1 



(3.29) 



Note that the total self-energy £ and vertex D are real and have two contributions of 
different origin: one arises from the interaction of the system and the bath (rj and v) and 
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one is caused by the non-linearities stemming from the average over disorder (that we called 
E and D in Eqs. ( |3.27| ) and ( 3.28Q ). If p = 2, X and D are identical to the classical ones. 
Instead, if p > 3, the non-linear terms acquire an explicit dependence upon h. 
The dynamic equations can then be written in a compact form 

/•CO 

(md 2 + z(t))R(t, t') = 5(t -t')+ / dt" m, t")R(t", t') , (3.30) 

Jo 

POO ft' 

{md 2 + z{t))C{t,t') = / dt"m,t")C(t",t')+ dt"D{t,t")R{t',t") , (3.31) 

J o Jo 

that we call later the i?-eq. and C-eq., respectively. It is important to realize that the 
self-energy Y,(t,t') is proportional to the response function R(t,t'), which in turns implies 

E(t, t') = E{t, t') = for t < t' . (3.32) 

This means that the upper limit of integration in Eqs. ( |3.3(J| ) and ( |3.31| ) is t, which renders 
the equations explicitly causal. There are no more independent equations for R and C . The 
other two equations that can be obtained from Eq. Q3.17|) are the equation for R(t',t), that 
is equivalent to the R-eq above, and one equation that identically cancels by virtue of the 
identity ( gig) . 

In their integrated form as Schwinger- Dyson equations the dynamic equations read: 

R(t, f) = G Q {t, t') + t dt" C dt'" G (t, t") £(t", t'") R(t'", t') , 
Jt' Jf 

C{t, f) = f dt" f dt'" R{t, t") D(t", t'") R(t', t'") , (3.33) 
Jo Jo 

with the propagator given by 

G~ x {t,t') = md 2 t + z(t) . (3.34) 

Real and imaginary parts of Eqs. ( 3.181 ) and (|3.19|) combined with Eq. ( [3.15 ) imply the 
equal-times conditions 

C(t,t) = l, 

R(t,t) = 0. (3.35) 

In addition, from Eq. ( |3.20| ) one obtains that the first derivative of the response function is 
discontinuous at equal times: 

lim d t R(t, t') = — , 
t'-*t~ m 

hm + d t R(t, t') = , (3.36) 

while from Eq. (|3.21|) one obtains that the correlation is continuous: 

lim dtC(t, t') = \im d t C(t, t')=0. (3.37) 
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The equation for z(t) = z + (t) = z~(t) can be derived from the Schwinger- Dyson equation 
( 3.33Q , by imposing the spherical constraint through the evaluation at t = t'. Multiplying 
operationally by G^ 1 one obtains 



z(t) 



dt" [E(t, t")C(t, t") + D(t, t")R(t, t")} 

+m f dt" f dt"'(d t R(t,t"))D(t",t"') (d t R(t,t"')) . 
Jo Jo 



(3.38) 



The last term is a consequence of having a kinetic term with second derivatives. Had we 
chosen a first derivative term in the kinetic energy the last term would have not appeared. 
It can be easily identified with minus the second-derivative of the correlation at equal times 
by taking the limit t' — > t~ in Eq. ( |3.31|) . Thus 



z (t) = f dt" Mt, t")C(t, t") + D(t, t")R(t, t")] - m ^C(t, t' 
Jo ot 2 



(3.39) 



In conclusion, Eqs. (|3.30|) , (|3.31|) and ( p.39|) are the complete set of equations that 
determines the dynamics of the system. In the following we make a time reparametrization 



that implies 



and 



C = C 



t = Mj t 



R 



M 7o 



R 



J=J 



f3 = /3 



A 



A 



M 7o 



m = (M7 D ) 2 m , U = M^ h , 

The system of units we use is 

[C] = [K] = [h] = i, 

[m] = 0} = [1/J] = [1/A] = [t\ . 

Hereafter we write the equation in reparametrized variables and drop all hats. 



(3.40) 



(3.41) 



(3.42) 



(3.43) 



IV. SOLUTIONS TO THE MODE COUPLING EQUATIONS 

As already mentioned, the set of integro-differential equations ( |3.30| ), ( |3.31| ) and ( |3.39| ) 
is causal - as in the classical case - and one can attempt a numerical solution by a simple 
algorithm. The correlation C(t,t') and response R(t,t') are two time quantities, that is, 
they depend on t (which physically corresponds to the time of observation) and t' (which 
corresponds to the age of the system). In a normal situation, one expects that after an 
equilibration time t EQ (in general model dependent) all two-time functions will depend only 
on the time difference t — t' . In other words, all two-time functions become time-translation 
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invariant (TTI). Under these circumstances, to which we refer as equilibrium evolution, 
response and correlation are not independent quantities but are linked via the fluctuation- 
dissipation theorem (FDT). 

Nevertheless, it is also known that certain systems never reach (at least for experimentally 
accessible time-scales) an equilibrium dynamics. For arbitrary large t' and in a certain region 
of the parameter space, although one-time quantities, such as the energy-density, tend to an 
asymptotic limit, two-time quantities, such as the correlation and response, depend both on 
t and t' in a non-trivial way. This is indeed the typical situation for a glassy system.BS 

The question we would like to explore is whether Eqs. ( |3.3U| ), ( |3.31| ) and Q3.39| ) encode 



a non-equilibrium evolution as described above. We know that indeed this is the case for 
classical disordered systems. As shown in Ref. |52[| , this kind of equations describe, at 
least partially, the non equilibrium dynamics of glassy systems, notably aging effects of 
spin-glasses. The salient feature of these solutions is that, below a critical temperature 
Td = T c (h = 0), and for large t', there exist at least two time scales: for t — t' small with 
respect to a characteristic time Tit'), that depends upon the age of the system, the dynamics 
is similar to an equilibrium evolution (stationary regime) but for t — t' large with respect to 
T(t') the dynamics becomes extremely slow with aging effects. In this regime, correlation 
and response are related in a manner that violates FDT. 

In this paper we show that this also happens in the quantum problem in a region of 
parameter space. In what follows we concentrate on the case p > 3. As in the classical 
case, the p = 2 model is connected to O(N) models and has a different behavior. In Fig. [1 
we present a schematic phase diagram (T, H) for p = 3. As we show below, the dynamics 
in the disordered ("paramagnetic') and ordered ("spin-glass') phases are characterized by 
"equilibrium' and "out-of-equilibrium" effects, respectively. The dynamics for the quantum 
system in the glassy phase has a similar pattern to the classical one, in the sense that 
there are at least two time-sectors in which the evolution is of a different nature. Inspired 
by the classical problem,0 we propose that also in the quantum case the weak-ergodicity 
breakingi'0 and the weak long-term memoryS properties hold. These two assumptions are 
based in part on the insight coming from the numerical solution of the full equations. 

A. The correlation function and the weak-ergodicity breaking scenario 

The weak-ergodicity breaking scenario states that, for t > t', the correlation function 
decays in such a way that 

lim C(t,t') = q + C ST (t-t') (4.1) 
lim C ST (t-t')=0 lim lim C(t,t') = q (4.2) 

t—t'—*oo t—t'^oot'—+oo 

lim C(t,t') = . (4.3) 

In these equations q is the Edwards- Anderson order parameter that characterizes the spin- 
glass phase. This means that for large t > t' such that t — f is small with respect to T[t') , 
the correlation function first decays from 1 to q in a time-translational invariant manner. 
T(t') is a growing function of t' whose precise form depends on the model. The correlation 
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PM 




FIG. 1. The schematic phase diagram for the p = 3 model. 

goes further below q to eventually reach in a manner that depends both upon t and t' (the 
aging effect). This behavior suggests the presence of at least two time-sectors in which the 
dynamics is stationary and non-stationary, respectively. 

We then write the correlation as the sum of a stationary and an aging contribution: 

C(t,t') = C ST (t-t') + C AG (t,t') . (4.4) 

The matching conditions at equal times between C ST and C AG are 

C(t,t) = l C ST (0) + C AG (t,t) = l (4.5) 

with 

C ST (0) = l-g C Aa (t,t) = q. (4.6) 

Together with Eq. (|4.2|) they ensure that in the two-time sector in which C ST decays from 
1 — q to 0, C AG is just a constant q. Instead, in the time-sector in which C AG decays from q 
to 0, C ST vanishes identically. 

The name weak-ergodicity breakingP0 reflects the fact that for short time-differences 
the system behaves as if it would be trapped in some region of phase space of "size" q - 
suggesting ergodicity breaking. However, it is always able to escape this region in a time- 
scale T{t') that depends upon its age t' . Hence, trapping is gradual and ergodicity breaking 
is weak. 
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We recall that in the classical and purely relaxational case (Langevin dynamics) the 
correlation functions are monotonic with respect to both times t and t'. However, in the 
presence of inertial terms there can be oscillations and the decay can be non- monotonic. 
This will depend upon the relative value of the mass m with respect to the other parameters 
in the problem. For the parameters we study in this paper, the oscillations appear only in 
the stationary regime, the aging dynamics having a monotonic decay towards zero. This is 
relevant since it will allow us to use the general properties of monotonic correlation functions 



proven in Ref. |12|1 to find the two-time scaling of C AG (t,t') 



B. The response function and the weak long-term memory scenario 

Regarding the response function, we propose 

R{t, = R ST {t - t') + R AG {t, (4.7) 

with 

R ST (t-f) = lim R(t,t') . (4.8) 

When a system is in equilibrium, the response is simply related to the correlation via FDT. 
Here we need to extend this relation to a non-equilibrium quantum situation. As in the 
classical case we assume that the dynamics in the stationary regime is like an equilibrium 
dynamics and satisfies FDT. We then assume that R st (t) and C ST (r) are related by 

R s Ju,) = -I hm / — tanh ( C ST (u) (4-9) 

(see Appendix ||). In contrast, in the time-sector in which the dynamics is non-stationary 
and manifestly non-equilibrium, FDT will be modified. In Section [VI]| we suggest an exten- 
sion of FDT to quantum non-equilibrium systems with slow dynamics, that we later check 
in the particular model under study. 

The weak long-term memory scenario statesEE that the system keeps a weak memory 
of what happened in the past. More precisely, the response function tends to zero when 
times get far apart, 

lim R(t, t*) = V fixed t* , (4.10) 

and its integral over a finite time-interval also vanishes 

but, its integral over an interval that grows with time gives a finite contribution: 



lim / dt"R(t,t") = V fixed t* (4.11) 

t— »oo Jo 



lim / dt" R(t,t") < Hm / dt" R(t,t") V fixed t* . (4.12) 

t-t* ^oo J t _ t * t^oo Jo 

Note that the property Q4.1(Jp implies 
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lim R ST (t-t') = 0. (4.13) 

t— V— >oo 

An important quantity that reflects the decay of the response function is the two-time 
susceptibility defined as the integral of the response over the time-interval [£„,,£]: 

X (t,t w ) = f f dsR(t,s) (4.14) 

J t w 

(/ is the strength of the applied perturbation, see Eq. (|2.8|)). In the classical case, a use- 
ful description of the departure from FDT is obtained by plotting x(Mw) vs C(t,t w ) for 
several t w using r = t — t w as a parameter. Since the classical FDT implies x(Mw) — 
l/T(C(t, t w ) — C(t, 0)), violations of the theorem manifest as departures from a straight 
line plot. In the quantum case, the relation between R and C is more involved and a priori 



one does not expect these plots to be useful. As we show in Section [VI \ it turns out that 
even for a quantum system these parametric plots are relevant. 



V. DYNAMICS IN THE PARAMAGNETIC PHASE 

In the paramagnetic phase (see Fig. [I]) one expects that, after a short non-equilibrium 
transient, a time-translational invariant (TTI) solution establishes. We show that a TTI 
ansatz implies that FDT holds between R and C and between E and D. By solving nu- 
merically the full set of equations (|3.3U|) , (|3.31|) and (|3.39|) we show the correctness of this 



ansatz and later discuss some properties of the solutions. 

After the short initial transient, one-time quantities such as z(t) should reach a limit z^. 
We evaluate it as 

/oo A^j 
— \E(u)C(u) + D(u)R(u)] -md 2 C(t-t') . (5.1) 

-oo 2ll t'—>t -^oo 

The Fourier-transform of Eqs. ( |3.30| ) and ( |3.31[ ) yields the quantum mode- coupling equations: 

R(lu) = = , (5.2) 

—mcu 2 + Zoo + 4:i](cu) — T,(u) 

C{uj) = (2hu(u) + D(u))\R(uj)\ 2 . (5.3) 



Thus, 



JmR(u) = -Im(S(cj) - An{uo))\R{uo)\ 2 , (5.4) 
ReR(u) = (-mcu 2 + z (X3 -Re - 4r}(uj))) \R(co)\ 2 . (5.5) 



If we now assume that E and D are related by FDT, it is easy to check that Eqs. Q5.2| ). 
( |5.3| ) and ( |5.4j ) imply that R and C are also related by FDT. The proof is completed by 
showing that FDT between R and C implies FDT between E and D. This is easier to do in 
the time-domain, by verifying that Eq. ( |B.8| ) holds between E and D if it does between R 
and C. 
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C(r + t w , t w ) 
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2 4 6 8 



r 

FIG. 2. The auto-correlation C(r + t w ,t w ) as a function of the time-difference r for 
= 2, 4, 8, 16 for the classical model in the paramagnetic phase, T = 3 and % = 0. The curves for 
tu, = 8 and ^ = 16 are superposed. As in all following figures, p = 3, A = 5, J = 1 and m = 1. 




2 4 6 8 

T 



FIG. 3. The same parameters as in Fig. ^. The response function R(t + t w ,t w ) as a function 
of the time-difference r for t w = 2, 4, 8, 16. TTI establishes for t w > 8. In the inset, the integrated 
response x( T + tw,t w ) vs. C(r + t w ,t w ) in a parametric plot for waiting time t w = 4 and r in 
[0, 28]. The classical FDT prediction is represented by the straight line of slope — 1/T = — 1/3; the 
X vs. C curve coincides with it showing that the classical FDT holds. 
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FIG. 4. The auto-correlation C(r + t w ,t w ) as a function of the time-difference r for 
t w = 2, 4, 8, 16 for the quantum model in the paramagnetic phase, T = and fi = 6. The curves 
for i TO = 8 and t w = 16 fall on top of each other demonstrating TTI. 

0.4 



0.2 



R(r + t w ,t w ) 



T 



-0.2 



~~ I 
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1 -0.5 0.5 1 

C{t + t w ,t w ) 



FIG. 5. The response R(r+t w , t w ) as a function of the time-difference r for the same parameters 
as in Fig. ||. The waiting times are t w = 2, 4, 8, 16 and TTI holds for t w > 8. In the inset the 
integrated response x( T + t w ,t w ) vs. the auto-correlation C(r + t w ,t w ) in a parametric plot for 
t w = 8. Since the system is in the paramagnetic phase with pure quantum fluctuations, the \ vs - 



C plot is not expected to give us further information (see, however, Section VII) 
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A solution of this type exists only in a restricted region of the phase diagram. In order 
to demonstrate the existence of a dynamic transition, we solve the full dynamic equations 
numerically and find that at a critical line, both the TTI and FDT assumptions break down. 

For the purpose of solving Eqs. Q3.3C L (|3.31|) and (|3.39|) numerically, we proceed as 
in the purely relaxational classical caseJ^J" Notice though that the presence of non local 
kernels rj and u, that resemble the delta-function and appear convoluted with the correlation 
and response, renders the numerical solution harder. The larger the cut-off A, the smaller 
the iteration step 5 we need to compute these integrals with a good precision. We found 
an acceptable data collapse for 5 = 0.0025, 0.005, 0.01 if A = 5. Note that the values of 
5 typically used in the purely relaxational case are at least one order of magnitude larger. 
This is the reason why the time-window we explore here is rather narrow though it suffices 
to show the trend of the solution. 

Even after rescaling time and all parameters in the problem, there are still many free 
parameters left: (J, p, m, A, h, T). We here concentrate on the dependence upon only a small 
fraction of this set. We fix hereafter p = 3, J = 1, m = 1 and A = 5 and study the behavior 
upon the remaining parameters T, h. We discuss the data for the correlation, response and 
integrated response for three points in the paramagnetic phase: 



• T = 3 > Td, h — 0. This corresponds to the classical problem with colored noise and 
inertia. The classical critical temperature for the purely relaxational case is ~ 0.6. 
In Fig. |2| we show the correlation C(r + t w ,t w ) vs. the time-difference r for t w = 
2, 4, 8, 16. In Fig. |3| we plot the response R{r + t w , t w ) vs r for the same waiting times. 
In both cases the decay is very fast, typical of upper critical dynamics, and TTI quickly 
sets in (cfr. the curves for t w = 8 and t w = 16 that are indistinguishable in both plots). 
In the inset, we plot the integrated response x( T + t w , t w ) vs. C{t + t w , t w ) for t w = 4 
using r as a parameter. The x vs - C plot very soon becomes a straight line of slope 
— 1/T = —1/3 showing that the classical FDT is satisfied. 

• T = 0, ft = 6 > h c . This corresponds to a strong quantum regime without thermal 
noise. The data for correlations and response in Figs. ^ and [| are qualitatively similar 
to the ones for the classical problem at high temperatures in the sense that an equilib- 
rium regime is quickly attained. The decay is fast and TTI holds. In the inset of Fig. |5| 
we display the x vs - C plot. In this strong quantum case at T = 0, the curve severely 
deviates from the classical expectation of a straight line of infinite slope. The quan- 
tum fluctuation-dissipation theorem is verified in Fig. |6| by comparing the response 
function obtained numerically with the convolution of the kernel tanh(/?7iu;/2) with 
the numerical correlation function, see Eq. ( |B.12|) . These two curves should be equal 
in equilibrium, as imposed by FDT, and we here check that this is indeed the case in 
the paramagnetic phase. 

• T = 2,ft=l> h c (T — 2). This is the quantum problem with both thermal and 
quantum fluctuations in the paramagnetic phase. In Fig. [7] we check FDT for these 
parameters and observe, in the x vs - C curve displayed in the inset, that the FDT 
becomes the classical. 



25 



0.4 



0.2 

R{t,t-r) 





-0.2 

2 4 6 8 



FIG. 6. Check of FDT for T = and fi = 6 (the paramagnetic phase). We compare the response 
R(t, t — t) for total time fixed t = 32 and r 6 [0, i] obtained from the iteration of the dynamic 
equations, and the response computed from Eq. ( p.l2j ) that assumes FDT. The correlation values 
used in Eq. (pUp are also obtained from the numerical algorithm. 

We continue the study of the paramagnetic phase by analyzing the variation with T and 
h, close to the transition line, of the out of phase susceptibility 

X"H = . (5.6) 

This is the quantity that is most commonly measured experimentally. For instance, in the 
context of the study of quantum phase transitions, the dynamics of the randomly diluted, 
dipolar coupled, Ising magnet LiHo x Yi_ :r F4 was studied by measuring,il among other quan- 
tities, 

We follow two routes to the transition: 

• T = 1.3 > Td, h = 3, 2, 1, 0. One approaches the classical model by decreasing the 
strength of quantum fluctuation at fixed temperature T — 1.3 ~ In Fig. || we 
plot }("{u)) vs. uj. In the curves for small h we observe a local minimum that is 
associated to the appearance of a plateau in the correlation function, at the value 
q, when approaching T^. This is typical of "discontinuous transitions" as the one 
in the classical p-spin modeHlJ (and in super-cooled-liquids B). We here note that 
this feature survives the presence of small quantum fluctuations and progressively 
disappears when h increases. For h = 3, the minimum has been substituted by a 
single peak. 

• T = 0.1 < T d , h = 6,5.5,5,4.5,4,3.5. One approaches the critical line (T c ,h c ) 
from above by modifying the strength of the quantum fluctuations while keeping the 
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FIG. 7. The same comparison as in Fig. ^ for T = 2 and fi = 1 (the paramagnetic phase). In 
the inset, the x vs C plot for t w = 16. The straight line of slope —1/2 is the expectation for a 
classical FDT at this temperature. For t w larger than a characteristic time, the FDT becomes the 
classical one, as expected in a high-temperature situation. 

temperature fixed T = 0.1. Figure |9| displays yf'{oj) vs. to for this series of h. (We 
cannot get closer to the transition line because the dynamics gets so slow that the 
stationary limit cannot be reached.) The curves have a distinct peak that moves 
towards decreasing frequencies for decreasing %. The spectral width of the bell-shaped 
curves also increases when approaching the critical line. The low frequency part of 
the curves is importantly suppressed with increasing fi. All these features show that 
quantum fluctuations radically affect the relaxation. 

Contrary to Fig. [8], there is no local minimum in any of the curves. This is in accordance 
with a decreasing and eventually vanishing plateau value q when the parameters ap- 
proach the quantum critical point. We shall further discuss the behavior of the model 



close to the transition line in Section VII B 



These results demonstrate the existence of a cross-over in the disordered phase, separating 
a region where the scaling laws are controlled by the quantum critical point (0, h c ) from 
another region with different scaling laws controlled by the classical critical point (T^,0). 

Before concluding our analysis of the paramagnetic phase, let us compare our results 
with the situation encountered for the system!! LiHo x Yi_ x F4, which is believed to be an 
experimental realization of the transverse Ising model.0 Naturally, we do not expect to find 



Note however that experiments signal a rather different dependence of the dynamics on the 
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the same behavior in both models as we know that in the classical limit they are different :E3 
the classical transition being first order for the p-spin model and second order for the Ising 
case. Our analysis indicates that the dynamic transition of the p-spin model remains first 
order along the critical line but becomes second order at the quantum critical point. This is 
the mirror situation to what has been observed in the dipolar glass. In Ref. p3j|, it has been 
suggested that the quantum spin-glass transition may be first order. Figure |8| displays the 
effect of quantum fluctuations close to the classical critical point and looks rather similar 
to Fig. 2 in Ref. [ 23 ] which instead shows the dynamics close to the quantum critical point. 
Figure || looks similar to Fig. 4 in Ref. [2^] , where the dynamics close to the classical critical 
point is presented. This comparison indicates that this model belongs to a different class 
from the system LiHo x Y 1 _ I .F4 and probably from the transverse field Ising model. 



VI. DYNAMICS IN THE GLASSY PHASE 

In this section we focus on the non-equilibrium dynamics observed in the glassy phase. In 
order to justify the scenario discussed in Section |T7], we first present the results from the nu- 
merical solution of the full equations. Next, we analyze these equations in the stationary and 
aging regimes and obtain a self-consistent relation that determines the Edwards-Anderson 
parameter. 



A. Numerical results 

The numerical solution shows a very different behavior below the critical line h c (T c ). In 
order to illustrate the slow dynamics, the breakdown of TTI and the weak-memory of the 
system we choose the following parameters: 



T < Td, h = 0. This is the classical, glassy model studied in Ref. JTT[| but with the 
addition of inertia. The effect of inertia amounts to the presence of oscillations around 
q but it does not change the qualitative behavior of the model. For example, as we 
shown analytically below, the value of q is not modified by inertia if m is small. 



T = 0, h = 1. This corresponds to pure quantum fluctuations. Figure [H] shows the 
correlation C{j + t w , t w ) vs. r for the waiting times t w = 0.5, 1, 2, 4, 8, 16. The curves 
have a stationary and an aging regime separated by the Edwards- Anderson parameter 
q ~ 0.72, see Eq. (f4.2|) . The correlations decay to zero at far apart times, as in 
Eq. ( |4 . 3| ) . This demonstrates the weak ergodicity breaking scenario. 

In Fig. [LI] we plot the response R(t + t w ,t w ) vs. r for the same waiting times. One 
observes a small time-difference regime, r < 5, where TTI rapidly establishes and a 



temperature history of the sample in the case of classical orientational glasses and classical spin- 
glasses with short-range interactions .i@ This might give a warning on choosing the transverse Ising 
model to describe the dynamics of quantum dipolar glasses. 
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FIG. 8. x"M vs. lo for T = 1.3 and ft = (O), ft = 1 (+), ft = 2 (□), ft = 3 (x] 
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FIG. 9. x"(w) vs - ^ for T = 0.1 and fi = 3.5 (★), h = 4 (A), ft = 4.5 (O), ft = 5 (+), ft = 5.5 
(□), ft = 6 (x). 
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large time-difference regime, r > 5, where the response shows a decaying queue with a 
t w dependence. The integration of this queue over a time-interval scaling with t w will 
tell us about the violations of FDT in the quantum case (see Fig. [16] below). It also 
shows that the system has a weak long-term memory as described in Section [TV| . 

• T = 0.5, h = 1. The qualitative behavior is similar to that obtained for T = and 
h = 1 apart from the fact that thermal fluctuations slightly decrease the value of 
the Edwards-Anderson parameter and accelerate significantly the decay in the aging 
regime, as demonstrated in Fig. 0. 



B. The Lagrange multiplier z(t): a one-time quantity that reaches a limit 

The Lagrange multiplier z(t) is a one-time quantity, i.e. it depends only upon the total 
time t. In this formalism we assume that all one-time quantities reach a well-defined limit 
asymptotically. In Appendix [B| we describe how we obtain the following equation for the 
asymptotic limit of z(t): 

roo roc 

z 0O = A o0 + q dr' S ST (r') + D q / dr'i? ST (r') 
Jo Jo 

+ / dr' [ £ st (t')C st (t') + D ST (r')R ST (r') ] - m d 2 T C ST (r) . (6.1) 
Jo 

Let us discuss each term in this expression. We called the aging contribution: 

Aoo = lim f dt" \ S AG (t, t")C AG (t, t") + D AG (t, t")R AG (t, t") 1 . (6.2) 

t^co Jo I J 

The kernels v and rj, that are related to the bath, do not contribute to A^. This is so 
because we have assumed that for r > T(t') the kernels have already decayed to zero in 
such a way that they do not contribute to these integrals. More precisely, we are neglecting 
terms of the form 

lim /* dt" A(t-t")B(t,t") (6.3) 

t^oo Jq 

where A is either v or rj and B is either C AG or R AG . 

The second and third terms come from the constant (non-zero) limit of the first decay 
of the correlation q = limi_i/^oo linn/^oo C(t, t') and equivalently of the vertex 



D q = lim lim D(t,t') . (6.4) 

In the model under study 

D q = ^ <f - 1 (6.5) 
if we use lim T ^ 00 -R st (t) C g, a property of the weak long-term memory scenario. 
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FIG. 10. C(t + t w , t w ) vs. r for T = 0, % = 1 and from bottom to top t w = 2.5, 5, 10, 20, 30, 40. 
The weak ergodicity breaking scenario and aging are explicit. 
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FIG. 11. R(t + t w ,t w ) vs. r for T = 0, h = 1 and from top to bottom t w = 5, 10,20,30,40. 
The weak long-term memory scenario is explicit. 
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FIG. 12. C(t + t w ,t w ) vs. r for T = 0.5, h = 1 and from bottom to top ^ = 4, 8, 16, 32. The 
Edwards- Anderson parameter is slightly smaller and the decay in the aging regime is faster than 



in Fig. 10 



C. Dynamical equations in the stationary regime 

If (£,£') are such that C(t,t') > q, the discussion in Section [TV] implies C(t,t') = q + 
C ST (t — t') and R(t — t') = R ST (t — t'). The Schwinger-Dyson equation for R in this time 
sector reads 



(md 2 T + Zoo ) R st (t) = 6(t) + f 



dr'E ST fr - r'l-RsTfr'l 



(6.6) 



and it keeps the same form as in the high-temperature regime, apart from the fact that the 
constant has contributions from the aging regime. 
The Schwinger-Dyson equation for C reads 



(mdl + Zoo ) (q + C ST (r)) = A 00 + qJ^ dr' E ST (r') + D q j dr' R ST (r') 

/oo 
dr' [E ST (r + t')C st (t') + D ST (r + r')i? ST (r / )] 
-oo 

One can now Fourier-transform both equations 



(6.7) 



R*t(u>) 



1 



-mui 2 + z c 



E ST M 



(-muj 2 + z^ C st (uj) + Zoo q8(u) = (A^ + gS ST (u;) + D q R ST (uj)^ 8(uj) 



(6.8) 
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The formal solution to the equation for C ST is 

C ST (u) = (-Zoog + + g£ ST (^) + D q R ST (u)j 5(^)i? ST (^) + D st (uj) \R st (uj) \ 2 . (6.9) 

The first term on the right-hand-side has an imaginary and a real part. The imaginary part 
vanishes identically since, due to FDT, both lm.R ST (u) and ImS ST (w) are proportional to 
tanh (flfkj/2) which is zero at uj = 0. Concerning the real part of this first term, as we have 
assumed that C ST (r) goes to zero for r — > oo, we need to impose the self-consistent condition 

- Zoo q + + gS ST (cu = 0) + D q R ST (uj = 0) = . (6.10) 

This is the condition that fixes the Edwards- Anderson parameter. We shall find it again in 
the next section as the matching condition between the stationary and aging regimes. 
The final equation for C ST (u) is 

Cst(^) = -Dst(^)|-Rst(^)| 2 . (6.11) 

One can check that these calculations are consistent with the results from z^. Actually, 
the integrals in z-eq. involving the stationary parts can be evaluated with the help of the 

, and yield once again Eq. (pTTD). 
Similarly to the high-temperature case one can now show that FDT for S ST and D ST 
implies FDT for R ST and C ST . The remainder of the proof, i.e. to show that FDT between 
R ST and C ST implies FDT between E ST and D ST depends only upon the from of S ST and D ST 
as functions of R ST and C ST and is not modified from the one discussed in Section [V|. 



equations for R ST and C ST , Eqs. (p.8|) and (6.11 



D. Dynamical equations in the aging regime 

If we now choose the times t, t! to be well-separated so as to have C(t, t') = C AG (t, t') < q 
and R(t, t') = -Rag (t > t'), the weak-ergodicity breaking and weak long-term memory hypothe- 
ses allow us to throw the second time derivatives on the left-hand-side. We assume that their 
contribution is much weaker than the one of each of the integral terms on the right-hand- 
side. This is an assumption that we have to verify at the end of the calculation, once the 
solution for C AG and R AG is known. It corresponds to the over-damped limit: 

mdtC AG < Terms in the RHS c (t,t') , 

mdlR AG < Terms in the RHS R (t,t') . (6.12) 
Using the approximation described in Appendix 0, the R-eq in the aging regime becomes 

z^R AG {t, = S AG (t, *') / dr'R ST {r') + R AG (t, t') / dr' S ST (r') 

JO JO 

+ / dt"t AG (t,t")R AG (t",t') (6.13) 
Jt' 

and we call it the i? AG -eq. Similarly, the C-eq becomes 

/■OO /"OO 

ZooC AG (t } t') = C AQ (t, t') / rfr'S ST (r') + D AG (t, t') / dr'R ST (r') 

Jo Jo 

+ f dt" S AG (t, t")C AG (t", + f dt" D AG (t, t")R AG (t', t") (6.14) 

JO JO 

and we call it the C AG -eq. In all integrals we approximated E AG (t, t') ~ T, AG (t,t') and 
D AG (t,t') ~ D AG (t,t') since T)(t — t') and v(t — t') are very rapidly decreasing functions. 
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E. The Edwards- Anderson parameter 



The Edwards-Anderson parameter q is determined self-consistently from the matching 
of lim^oo C AG (t, t) = lim^/^oo linu/^oo C(t, t') = q. Taking the limit t' — > t~ in the i? AG -eq 
and C AG -eq one obtains 

_ roc roo 

ZooRAo(t, t) = £ AG (t, t) / dr' R ST (r') + R AG it, t) / dr' S ST (r') , (6.15) 

Jo Jo 

POO _ rOO 

z oc q = A OQ + q dr' S ST (r') + D AG (t, t) / dr' i? ST (r') . (6.16) 
Jo Jo 

The first equation admits the solution R AG (t,t) = since £ AG (£,t) is proportional to 
-Rag(M) ~ see Eq. ( p.27| ). This corresponds to the high-temperature solution where there is 
no aging regime. Here we concentrate on the other possibility, that is to say when 

S AG (t, t) f°° , 



-R AG (t, t) 



roc poo 

/ dr' R ST (r')+ / dr'E ST (T') . (6.17) 
Jo Jo 



In actual fact, we can further approximate this equation by using one of the assumptions 
already used all over this section: that the response becomes smaller and smaller as time 
passes - though its integral over an infinite interval gives a finite contribution. If we neglect 
all terms that are proportional to R AG (t, t) with respect to terms that are proportional to q, 
only the first term in the power expansions of £ and D survive: 

t/R) = l im ^M d = limD AG (M) (6.18) 

/q t^oo R AG {t,t) 

that in our model become 

t/R) = J2piP ~ 1] g p - 2 A^^V' 1 , (6.19) 



in accord with the large r limit of the stationary values ( |6.5| ). Equations ( |6.15| ) and ( |6.16|) 
become 

p oo p od 

Zoo = (t/R) / dr' R ST (r') + / dr' E ST (r') , (6.20) 
v J iJo Jo 

rOO _ poo 

z 00 q = A 00 + q dr' S ST (r') + D q / dr' R ST (r') . (6.21) 



The second equation is the same as the one arising from the end of the stationary regime, 
Eq. ( 16301) . 

We now study these equations with the help of Eqs. (|6.8|) and (|6.11| ) for the stationary 
functions. Using 

r°° l 

/ drR ST (r) = R ST (uj = 0) = —- -, (6.22) 

Jo Zoo - S ST [UJ = 



one has 
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FIG. 13. The ft dependence of the Edwards- Anderson parameter q that is read from the plateau 
in the log-log plot of the auto-correlation C(r + t w ,t w ) vs. t for t w = 32. From top to bottom 
ft = 0.1, 0.5, 1, 2 and the temperature is always T = 0. 



t/R] RUu = 0) . (6.23) 
i 




We remind that the factor Rs T (uj = 0) can be written in terms of the stationary correlation 
function using FDT; therefore this is a closed equation for the correlation that gives the 
Edwards-Anderson parameter q. In the case of the p-spin model it reads 

~^tanh(^W)Y. (6.24) 

-oo w' \ 2 J J 

In the classical case, the integral can be readily computed and the final equation for q is 

J2P{P 2 ~ 1) q^\l-qf=T\ (6.25) 

that coincides with the result for the purely relaxational dynamics.0 For p > 3 fixed, q is a 
function of temperature; it equals one at T = and tends to qd = q(Td) > at the dynamic 
critical temperature T d = T c {% = 0). 

If ft 7^ 0, the Edwards- Anderson parameter q depends upon ft and T, q(T, ft). It decreases 
when thermal and/or quantum fluctuations increase. One observes ^(0,0) = 1. At the 
critical line (ft c , T c ), q(T c , ft c ) ^ as long as T ^ 0. At the quantum critical point the nature 
of the transition changes and q(h c , 0) = 0. We shall explain this result after introducing the 



modification of the quantum FDT needed to complete the solution. Figure [13] shows the 
tendency of q(T = 0, ft) for ft = 0.1, 0.5, 1.2. 



35 



VII. FLUCTUATION-DISSIPATION RELATIONS OUT OF EQUILIBRIUM 



One of the important characteristics of the dynamic evolution in equilibrium is the 
fluctuation-dissipation theorem theorem (FDT) that establishes a model independent re- 
lationship (in the sense that only involves T and h but no other parameter in the model) 
between correlation and response functions. As recalled in Appendix |B[ this result can be 
easily obtained, for the quantum and classical cases, under the hypothesis of equilibrium. 
In general, such a relation is not expected to hold in a non equilibrium situation and we 
can then say that FDT is "violated" (of course what is violated is not the theorem but the 
hypothesis under which the result is derived) . It is nevertheless surprising that in the glassy 
phase of certain classical models, the violation of FDT0 takes such a simple form that it 
is possible to formulate a generalized version of this relation through the introduction of an 
effective temperature T EFF .E3 In this way the FDT takes the standard equilibrium form but 
with a model dependent effective temperature which is different from the bath temperature 
T (in general T EFF is higher or equal than T and it can take different values in different 
time-sectors). We would like to explore in this section how these ideas carry over to the 
quantum case. 



A. General considerations 

The standard quantum fluctuation-dissipation theorem states 

R(t - if) = % 9(t - t') [°° — exp(-iu(t - if)) tanh ( ^ j C(u) (7.1) 
n J-oo 2ix y 2 J 

with 

poo 

C(u) = 2Re drexp(iuT)C{T) . (7.2) 
Jo 

Guided by the classical case, we introduce the effective temperatureil 

T -^ s xm < 7 - 3) 

and we propose, as a generalization of FDT to the out of equilibrium case, the following 
relation 

R(t, t>) = ^6(t- t') r — exp(-tuj(t - t')) tanh ( X ^ t '^ hu ) C u u ) (7.4) 
n J-00 2n \ 2 J 

with 

C(t, u) = 2Re / ds exp(iw(t - s))C(t, s) (7.5) 
Jo 

and X(t,t') a function of both times t and if. Evidently, Eqs. (|7.4j) and (|7.5|) reduce to 
Eqs. ( |7. lj ) and ( |7.2| ) when i — > 00 if the evolution is TTI and the factor X is set to one. 
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Let us assume, for reasons that are intimately related to the time-separation discussed 
in Section [V, that X(t,t') is characterized by 



, fl if t-f<T(t') =► C(t,f) > q 

l ' > \X AG (T,h) if t-f>T(t') C(t,t') < q { ' 

For times in the stationary regime t — t' < T(t') this implies: 



R*r(t -f) = T 0(t - f° 7^ exp(-iu(t - t')) tanh ^ C ST (u) (7.7) 



with 



C ST {u) = 2Re / dr exp{iujT)C^{T) . (7.8) 
J o 

If we suppose that (even at T = 0) /3X AG (T, h) is finite in the glassy phase, it is then clear 
that in the aging sector the Fourier integral is dominated by u ~ 0. After using 

tanh (*f£\ . jULjOg!! (7 . 9) 

we obtain 

R AG (t, t') = (3X AG (T, h) ^-C AG (t, t') . (7.10) 

at 

For a standard quantum system in equilibrium at finite temperature, there is a char- 
acteristic time, that depends upon temperature and the strength of quantum fluctuations, 
beyond which the evolution becomes classical. Concerning FDT, this means that the clas- 
sical relation establishes with X(t,t') = 1. This is clear, for instance, for the model under 
consideration at finite temperature in the paramagnetic phase (see the inset in Fig. ^). If 
our assumptions are correct, the situation in the glassy phase is different. For small times 
compared with T(t w ), now a waiting-time dependent characteristic time, correlation and 
response are related by the standard (quantum) FDT. Instead, for time differences larger 
than T(t w ), they are related in the same way as in the classical case but with a T and h 
dependent effective temperature that is different from the one of the environment. 

In general, the proposal given in Eq. ( |7.10| ) will allow us to reduce the C AG -eq. and R AG -eq 
equations to single equation for C AG . It is simple to check that, at least for the model under 
consideration, this equation has the same structure as the one for the classical problem. All 
dependence on h enters through q, X AG and R ST . Hence, once the explicit form of X AG is 
known, the aging equations are solved by the same ansatz as in the classical case: 

C*M = V- 1 (t&] , (7-11) 




R AG (t,t') = ^qd t/J - ] j^-l . (7.12) 
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Note that the general arguments introduced in Ref. ||12|| to characterize the behavior of 
monotonia correlation functions apply to the correlations in the aging regime and justify the 
solution (|7.11|) - ([7.12 ). Even if the inertia might make the correlation oscillate, it does only 
in the stationary regime since it oscillates around q, and it becomes monotonic in the aging 
In Eq. (|7.6|) , we assumed that X(t,t') can take only two values in view of the model we 
study However, in the classical case, a more general situation arises in models with multiple 
time-scales, for which the factor measuring the violation of FDT is a non-constant function 
of the correlation in the aging regime.0'u We expect the same structure to carry through 
to quantum problems whose classical counterparts are of the multi-scale type. 



B. Analysis of the p > 3 model 

We check, both analytically and numerically, the viability of our assumptions for the 
p-spin model. The factor X AG measuring the violation of FDT is determined by Eqs. (|6.21 ) 



and ( ET22D 



= A, 



+ D q R ST {u = 0) . 



R st (lo = 0) 

Using Eq. ( [7.1 0|) and the equivalent relation between E AG and D AG , we obtain 



(7.13) 



and 



^ lim (D AG (t,t)C AG (t,t)) = ^fqD q 



T t^oo 



(7.14) 



(P-2) 



T 



R st (lu = 0) 



2(P - 2) 2 g - P /2 



\ Pip ~ 1) 



(7.15) 



In the classical limit X AG = (p — 2)(1 — q)/q and the result in Ref. [[TT| ] is recovered. Note 
that both in the classical and quantum case, X = if p = 2. Since the case p = 2 is formally 
connected to ferromagnetic domain growth (in the mean-field approximation) there is still 
no-memory in the quantum domain growth. 

The classical critical point (T^, h — 0) is characterized by§ 

X AG (T d ,0) = l g(T d ,0)^0. (7.16) 

Following the critical line towards the quantum critical point (T = 0, h c ), one obtains 

/ 2(p- 2? N 1/p 



-^Ag{T C j foe 



q(T c , H c 



\J 2 p(p — 1), 



rp2/p 



(7.17) 



and a finite zero-frequency stationary susceptibility R st (uj = 0) < +oo. At the quantum 
critical point, q vanishes and the zero-frequency stationary susceptibility diverges: 



R ST {u = 0) ocT c ( 2 - p )/ p -> oo . 



(7.18) 
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R(t,t- t) 




FIG. 14. Check of the quantum FDT for the quantum model in the glassy phase when times 
take values in the stationary regime C > q. The parameters are T = and H = 1. Thin curve: 
R(t, t — t) from the numerical algorithm with fixed t = 32; bold curve: R(t, t — r) from Eqs. (7.7) 
and (|7.8| ) using C st (t) = C(t,t — t) — q with C(t,t — r) from the numerical algorithm and q ~ 0.7. 
The curves coincide for r < 7. 



(It is interesting to notice that the case p = 2 is clearly different.) On the other hand, if 
approaching the quantum critical point from below along the axis (T = 0, h) 



T EFF (T = o,h) ~ (h c -hy 



then 



g ~ (h c - %) 



2a /p 



o) ~ % - h) 



a(l-p/2) 



(7.19) 



(7.20) 



We next check numerically the relation between response and correlation in the glassy 
phase, T = 0, h = 0.1 < h c corresponding to pure quantum fluctuations. 

First, we use Fourier analysis to check the quantum FDT for short time differences. In 
Fig. [H] we compare two different ways of obtaining R(t, t — r) for total time t fixed and 
equal to 32. The thin curve is the result from the direct numerical solution of the dynamic 
equations and shows a fast evolution for short time differences and a tail that gives rise to 
the weak memory of the system. The bold curve is obtained from Eqs. ( |7.7| ) and ( |7.8|) using 
C ST (t, t — r) = C(t,t — t) — q with C(t,t — r) obtained numerically. We have estimated the 
Edwards- Anderson parameter as q ~ 0.7. One observes that for small time-differences the 
two curves coincide showing that FDT holds in this time-sector. 

Second, we check the generalized FDT in the aging regime. For this purpose, it is 



convenient to use the x vs - C plots defined in Section [iy|. In Fig. [15] we observe, as 
predicted by Eqs. (|7.7|), ( |7.S| ) and ( [7.10[ ), two behaviors according to the relative value of 
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0.6 0.8 1 

C(r + t w , t w ) 

FIG. 15. The x vs C plot for T = and ft, = 0.1. From bottom to top, different curves are 
associated to t w = 5, 10, 20, 30. One observes how a constant asymptotic limit is approached for 
increasing t w . The dots are a straight line of slope — 1/T EFF = —X AG /T = —0.6 (the analytic 
result). 

1 

0.8 

x(t + t w ,t w ) o.6 
0.4 
0.2 


0.2 0.4 0.6 0.8 1 

C(Y + t w , t w ) 

FIG. 16. The x vs C plot for T = and Ti = 1. The waiting times are t w = 5, 10, 20, 30 and 
the dotted straight line has a slope — 1/T EFF = —X AG /T = —0.9 (the analytic result). 
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0.2 0.4 0.6 

C (t -h t W jt Wj 



0.8 



FIG. 17. The x vs - C curves for the quantum model at /i = 0.1 and different temperatures. 
From top to bottom T = 0.1, 0.2, 0.3, 0.4, 0.5, 0.8. In the inset, the estimated effective temperature 
T EFF = T/X AG as a function of temperature for the same h. See the text for a discussion. 



C and q. The portion of the graph that is significant for the aging regime corresponds to 
C < q. As claimed in Eq. (|7.10|) the x vs C curve approaches a time-independent asymptotic 
limit that is a straight line of slope —X AG /T. In this case, T EFF = 1.6 and this shows that 
even a T = 0, T EFF ^ 0. 

The effective temperature depends on T and % in a non trivial way. In Figure 16 we 
show the x vs. C pl°t f° r ft =1 and T = 0. In this case T EFF = 1.1. 

In Fig. [17] we work at fixed ft = 0.1 and study the variation of T EFF with T. From the 
X vs C plots we determine how T EFF varies with T by measuring the slopes of the different 



curves and by using Eq. ([7.15 ). The six curves above the dotted straight line correspond 
to temperatures below T c . From top to bottom T = 0.1,0.2,0.3,0.4,0.5. One notices that 
T EFF > T in all these curves. The curve below the dotted straight line corresponds to 
T = 0.8 > T c . It has a slope — 1/T (indicated by the dots) and T EFF = T. In the inset we 
plot the temperature dependence of the effective temperature. T EFF is linear in T when the 
bath temperature is above T c but severely deviates from the linear behavior for temperatures 
below the critical line. Two curious features are that, in the glassy phase, T EFF increases 
when either T or h decrease. The surprising increase of T EFF as T decreases is also obtained 
in the classical limit .0 
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VIII. SUMMARY AND CONCLUSIONS 



In this paper we presented a formalism that allows us to analytically study the behavior 
of quantum glassy models within the typical situation encountered in experiments, i.e. out 
of equilibrium. 

Faced to the question as to how relevant quantum fluctuations are for the asymptotic 
real-time dynamics, one would have naively answered that they are irrelevant since quantum 
mechanics is not expected to play a relevant role in the large time sector, meaning large t w 
and for all t >t w . We show in this paper that this would have been a wrong conclusion and 
that the dynamic behavior is much richer: 

Quantum effects, together with temperature, determine where in parameter space a glassy 
phase exists. We show that sufficiently strong quantum fluctuations destroy the glassy phase 
at arbitrary low temperatures driving the system towards a paramagnetic phase. We predict 
then the existence of a dynamic quantum phase transition. This has no analogue in the 
classical model where glassiness is always present at T = 0. 

Quantum fluctuations also dictate when, in the glassy phase, non-equilibrium effects 
manifest. For any waiting-time t w after preparation, there exists a characteristic time T(t w ) 
that determines the end of a stationary regime and the entrance into the aging regime. This 
means that for time differences t — t w < T(t w ) the dynamics is stationary while for t — t w > 
T(t w ) two-point functions have explicit waiting-time dependences. The characteristic time 
T(t w ) depends on temperature, the strength of quantum fluctuations and the waiting-time 

tyj . 

The dynamics in the paramagnetic phase satisfies the fluctuation-dissipation theorem. In 
the glassy phase, we found that the quantum FDT is satisfied in the stationary regime and 
that a generalization is needed in the aging regime. This generalization defines a two-time 
function X(t, t') that is consistent with the definition of an effective temperature T EFF (t, t') = 
T/X(t,t') and that we later characterized, based on the numerical evidence for this model. 
When T EFF is different from zero, the modification of the quantum FDT takes a similar form 
as the classical one but with a T and h dependent effective temperature. 

The characteristic time T(t w ) is reminiscent of a decoherence time in the sense that it 
separates the stationary regime in which the influence of quantum effects is very strong 
and explicit from the aging regime in which the influence of quantum effects is implicit, 
appearing through, e.g., the h dependence of T EFF . 

One of the virtues of the approach we have followed is that it is close in spirit to the 
dynamic approach to classical systems. This encourages us to attempt to extend many of the 
recent advances in the understanding of classical glassy dynamics to the quantum problem. 
We mention below some of the questions that merit further attention. 

The model here studied has a modification of the quantum FDT given by a piecewise 
function X = 1 if C > q (T(t w ) > t — t w ) and a constant X AG (T, %) < 1 if C < q 
(T(t w ) < t — t w ). It is known that in classical systems a more general situation can arise. 
Another family of models exists with the factor measuring the FDT violation given by 
X = 1 if C > q and a non-constant function of the auto- correlation X AG (C) < 1 if C < q. 
We expect the arguments put forward in Section |V11| to apply to this kind of system with 
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minor modifications. It would be of particular interest to see how this structure appears in a 
concrete model as, for instance, the problem of a particle in a long-range correlated random 
potential. 

In the classical case, an intriguing connection between the factor that measures the FDT 
violationilHHilil in the non-equilibrium approach X AG (C) and the Parisi function x(q), 
< q < c[ea of the replica analysis in equilibrium, has been noticed and checked in several 
occasions. The structure of both functions is the same and, on top of this, their explicit 
values coincide if an argument of "marginality" is used to determine x(q) in the equilibrium 
case. In Ref. PBfl , several quantum models related, but not equal, to the one considered here 
were studied with a Gaussian variational method and replica theory in the imaginary time 
Matsubara representation. In order to obtain a physical value for the conductivity of the 
model, Giamarchi and Le Doussal used the marginality condition to fix the replica breaking 
point parameter x = x(q), q < qEA- A careful comparison of the dynamic X AG with the 
static x is now in order. Nevertheless, this cannot be readily done by comparing the result 
in Ref. |26| with ours since the way in which both systems are coupled to the environment 
are not necessarily equivalent. Thence, if a formal relation between FDT- violations and 
the replica x(q) exists at the quantum level is an interesting question that deserves further 
study. 

The last remark we wish to make about the FDT violations is that in the classical case 
it has been related to the production of entropy during evolution.i3Il It would be desirable 
to see whether such a connection also holds at the quantum level. 

Concerning the formalism here developed, it can be extended in several directions. On 
the one hand, one would like to study models with a spatial structure such as the one asso- 
ciated to the motion of a manifold of internal dimension d in a random quenched potential. 
This quantum problem has many interesting applications and, in the case of an embedding 
space with a large dimensionality can be treated, for example, within a dynamic Gaussian 
variational approach. Elil On the other hand, it would be extremely interesting to study the 
opposite limiting case of low dimensions and analyze, for example, the real-time dynamics 
of the Ising chain in a transverse field. In addition, one could investigate the effect induced 
on the dynamics by different environments like considering sub or super Ohmic baths, non- 
linear couplings and even baths of a completely different nature as those generated by spin 
variables.^ 

The Martin-Siggia Rose generating functional for classical stochastic processes takes a 
particularly simple form if the time-space is enlarged to a superspace by introducing a couple 
of Grassmann coordinates, and the variables of the problem are encoded in a superfield.0 
This can also be done for systems with disorder .0 The dynamic equations are then written 
in terms of a single super-correlator that encodes the auto-correlation and response. The 
equilibrium theorems, viz. causality, invariance under time-translations and the fluctuation- 
dissipation theorem, follow from the symmetries of the MSR effective action written in 
this way.0 The supersymmetric formulation of the classical dynamics has been helpful to 
find a solution of mean-field like models with an internal dimension like the problem of the 
manifold in a random media but, most importantly, it has highlighted the formal connection 
between the out of equilibrium formalism and the replica approach to the statics. It would 
be very interesting to identify a supersymmetric structure in the quantum problem, probably 
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at the level of the Keldysh-Schwinger generating functional. This could provide a means 
to establish a relation between the quantum dynamic formalism here presented and the 
quantum replica approach. 



We believe that an analytical study of aging effects along the lines presented in this paper 
will be relevant in the understanding of realistic glassy systems. 

Indeed, even at very low temperatures, where quantum effects are important, a glass 
is a non-equilibrium system that evolves over many orders of magnitude in timel^l Aging 
effects, demonstrating this non-equilibrium evolution, have been observed with, for example, 
burning hole experiments in organic glasses (see Ref. |5^] and references therein). In the 
context of magnetic disordered systems with explicit quenched disorder like LiHo :I .Y 1 _ a .F4, 
a great effort has been devoted to the study of the quantum phase transition and the 
dynamics close to the quantum critical point. It would be interesting to perform dynamic 
measurements in the glassy phase of these compounds to search for aging effects and check 
if a scenario as the one here discussed with, for example, a two step relaxation of correlation 
functions and FDT violations is also present. 

Another area of research where the relevance of quantum effects in the very low tem- 
perature dynamics is currently under study are systems of magnetic nano-particles such as 
7— Fe2C>3 particles (maghemite) embedded in a silica matrix. The search for a non- stationary 
evolution in samples at small concentrations (~ 0.3%) has given a negative result.B However, 
in experiments performed at higher temperatures (where quantum effects are not expected 
to play an important role) with samples with a larger concentration (~ 17%) much of the 
phenomenology of spin-glasses has been recovered.E2l It is then licit to ask whether there is 
a region in the temperature-concentration parameter-space where the systems exhibit aging 
driven by quantum fluctuations. In the area of magnetic nano-particles, a notion of an effec- 



tive temperature greater than the temperature of the environment was advocated in Ref. [74 
and it would be interesting to explore the connection with effective temperature defined in 
this paper. 

We have chosen to comment on three kind of materials where effects like the ones obtained 
in this article might be observed. Surely enough, many other glassy systems with or without 
explicit disorder might exhibit a similar phenomenology. 
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APPENDIX A: THE SCHWINGER-KELDYSH PATH-INTEGRAL 



The generating functional ( |2.1| ) is more conveniently expressed with the help of path- 
integrals: 

Z[t,C] = JdcP (0|T*exp (-1 jHdf£-(f)0(t))Texp (~ jT dtt{t)<m) IflPM 
= J dWfadfa (0|T*exp (-y™dtC(t)<f>(t)) \<t>i) 



(A.l) 



The two matrix elements are now expressed as usual Feynman path-integrals 



d(pd<p 1 d<p 2 



"(OO) = 1 



0+(O)= 



"(oo)=i 
0"(O)=0 



P0" 



exp 



<0 2 |/3(O)|</>>. (A.2) 



The first integrals / d(p, etc. are standard while we denote / V<p + the path-integrals. The 
doubling of degrees of freedom is a consequence of having two matrix elements, 

one for each source Formally this expression is rewritten as 



Z[t,C]= J Vcf> + J Vcf>- exp [~ (s[<f> + ]-S[<f>-] + Jdtt<f> + - Jdtr<f>' 
x(0 + |p(O)|0->. 



(A.3) 



APPENDIX B: QUANTUM FLUCTUATION-DISSIPATION THEOREM 



In this appendix we recall the quantum fluctuation-dissipation theorem. Proofs and 
descriptions of this theorem can be found in several textbooks.Ellzl We express it in the 
time-domain in a form that we use in Section [V| to show that the dynamic equations are 
compatible with a TTI-FDT ansatz. We also write it in a mixed time-Fourier notation that 
gives us insight as to how extend it to the case of glassy non-equilibrium dynamics. 

If at time t' the system is characterized by a density functional p(t'), the two-time 
correlation functions read 

C AB (t,t') = (A(t)B(t')) = ^Tr [A{t)B{t')p{t')\ (B.l) 

where the time-dependent operators, in the Heisenberg representation, are defined as 

O(t)=exp(^)o(0)exp(-^) . (B.2) 
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The trace is defined in the usual way, Tr[«] = J2a(^a\ • \i>a), with {ip a } an orthonormal 
basis in Fock space. The normalization is given by Z = Tr [p(£')]. 

Since operators do not commute, the quantum two-time auto-correlation functions are 
not symmetric in times 

C AA (t,t') = (A(t)A(t')) ^ (A(t')A(t)) . (B.3) 
One thus defines the symmetrized and antisymmetrized correlation functions: 

C {A ,B}(t, = ~ (A{t)B(f) + B{t')A{t)) , 

C [A , B] (t,t>) = i (A(t)B(f) - B(t')A(t)) , (B.4) 

respectively. 

In linear response theory R AB (t,t') (see Eq. (|2.8| ) for the definition) and the correlation 
C[A,B](t,t') are related by the Kubo formula 

R AB (t, = ~e(t- ( B(t')} ) = jO{t- t') C [A)B] (t, t') . (B.5) 

If the system has reached equilibrium with a heat-bath at temperature T at time t', the 
density functional p(t') is just the Boltzmann factor exp(— (3H). It is then immediate to 
show that, in equilibrium, time-translation invariance holds 

C AB (t,t')=C AB (t-t') . (B.6) 

In addition, Eq. ( |B.1| ) with p = exp(~j3H) imply the KMS condition 

C AB (t, t') = C BA (t\ t + i/ffi) = C BA (-t - iph, -t') . (B.7) 

Using now the KMS properties and assuming, for definiteness, that t > it is easy to verify 
the following equation 

C{a,b}( t ) + jRab(t) = C {AB} (r*) - jRab(t*) , (B.8) 

where r = t + i(3Ti/2. This is a way to express FDT through an analytic continuation to 
complex times that we use in Section [V] to show that a TTI and FDT solution holds in the 
paramagnetic phase. 

In terms of the Fourier transformed C AB (u) defined by 

/OO (jijj 
— exp(-iuj(t-t'))C AB (uj) (B.9) 

the KMS relations read 

C AB {u) = exp(phu})C BA (-u) (B.10) 
and lead to the following relation between Fourier transforms of the correlation functions: 
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2C[ A ,b](u) = (1 - exp(-(3huj)) Cab{u) , 
2C {a ,b}{uj) = (1 + exp(-/3W)) C ab (uj) , 

C [a ,b){uj) = tanh (^H C{a,b}(u) • (B.ll) 

Back in Eq. ( p.5| ) this implies FDT 

RA B {t-t>)= l -e(t-t>) r — exp(-iu(t-t')) tanh(^) C {A)B} (u) ( B - 12 ) 

ft J—oo 71 \ Z I 

Using 

/ rft exp(iu;t) = lim = tt5(uj) +i — (B.13) 

Jo e^o+ u + ie v ' uo v ' 



one has 



Rab(co) = -i lim f° — L— tanh C Rs} (o;') (B.14) 

ft e->o+ J-oo n u — u' + ie \ 2 J 11 



from which we obtain the real and imaginary relations 
lmR AB (uj) = i tanh ^) C {A ,B}{u') , 

ReR AB (u) = -^f- — ^ tanh (^-) C {AtBy {u/) . (B.15) 

ft J—oo 7V LO — Co 1 \ Z J 

If /3ftcu/2 < 1, tanh(/?ftu;/2) ~ /?ftcu/2 and Eq. (gig ) becomes the classical FDT: 

1 dC A B(r) 

«ab(t) = -- — ^ — (B.16) 

with t = t — t'. 

In the zero temperature limit, and for A = B, FDT reads 

itfw) = -- lim / sign(u/) C(u') , 

v ; ft e ^ + J-oo 2tx uj - u' + ie V ; V ; 

R(t) = -6(t) / — sin(ut)C(u) . (B.17) 

ft JO 71 

A finite integrated response 

r*jKt)=fl( U =o)=ir— ^ < b ' i8 > 

jo n Jo it to' 

implies C(u> — 0) = 0. In the glassy model we study in the paper, below h c (T c ), the 
response and correlation turn out to be the sum of two contributions. One is a stationary 
part that satisfies FDT. Since one expects it to yield a finite contribution to the susceptibility 
x( T + t w , t w ), one must have: 



CO 



drR ST (r) Finite / drC ST (r) = 0. (B.19) 



Thus, the stationary part of the correlation function oscillates around zero or, in other 
words, the full correlation oscillates around q at zero temperature. This is observed in the 
numerical solution to the dynamic equations. 
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APPENDIX C: THE CLASSICAL LIMIT 



The classical limit can be checked at every stage of the calculation. We shall here show 
how one can recover (i) the Martin-Siggia-Rose (MSR) action from the Keldysh-Shwinger 
(KS) action Q3.2p ; (ii) the equations of motion, for the classical correlation and response 
from the equations of motion for their quantum partners. 

This Appendix is intended for readers already familiar with the MSR formalism. We 
refer to Ref. |36|] for a detailed description of the method. 

Let us identify, in the limit (3% — > 0, the quantum linear combinations on the left with 
the classical variables on the right 

cr + (t) — cr~(i) A/s ,,. . . % A , . 

— ^ — - - m , <r+(t) -> s (t) + -s(t) , 

— ^ ^->«(t), cr-(t) - «(t) - -s(t) . (C.l) 

In the classical limit h — *■ 0, and when the cut-off tends to infinity, the kernels rj and i/ 
become 

Arj(t - t') = 4M6(t - £') 7o 5'(t - t') , 
2hv(t - t') = AM lo k B T5[t - t') . (C.2) 

With these identifications, we establish a connection between the KS and the MSR actions. 
The friction and thermal terms become 

% /*oo ^ roo 

lim -S T = -2M 7o / dt is(t)—s(t) + 2k B TMj Q / dt (is(t)f . (C.3) 
ft->o n Jo dt Jo 

By calling T^ 1 = 2M7 Q we recover two terms of the MSR action. 
The kinetic and constraint terms reduce to 

i „ r°° / „o l 



% f°° ( ! 

lim -S = dt -is(t)md?s(t) - -is(t)(z + (t) + z~(t))s(t) 
n^o n Jo \ 2 



A roo . . 9 TV r°° / \ 

+ zh L dts{t) ~ s{t) + h~2 L dt ( z+{t) ~ ^ (t) ) • (C ' 4) 

By calling z + (t) — z~(t) = z Q (t) and z + {t) + z~(t) = z(t) the last two terms impose the 
spherical constraint and are usually included in the path-integral measure, while the second 
term is the one left in the MSR action: 

lim -S = - dt is(t)(mdt + z(t))a(t) . (C.5) 
ft^o a Jo 

Finally, the terms that depend upon disorder become 

lim ~ r dt V[ct + , J] - V[tr~, J] = - I" dt is ^IflZl . ( C .6) 
R-o hJo 1 1 Jo 6s ' 
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Putting everything together, in the large friction limit 1 = 2M r ) a 3> 1 such that 
one neglects the second derivative terms in the kinetic part of the action, one recovers the 
usual MSR action for the Langevin dynamics of the classical model in a thermal bath at 
temperature T: 



lim — Sks = — dt 

h-+o n J 



is{t) (r- l j t + z(t)) s(t) - ksTT- 1 ^s(t)) 2 + ^^W) 



(C.7) 



In order to study the classical limit of the dynamic equations (|3.30|) and ( |3.31| ) we 
recognize the classical limit of the correlation and response as follows: 

C(t,t>) = \(*+{t)*-{t') + cr-{t)cr + {t')) -> (s(t)s(t')) , 

R(t,i/) = % -(<j+{t)a-{t') + a + {t)a-(t')) -> (s{t)is{t')) . (C.8) 
n 

In a similar way, it is now easy to verify that the classical dynamic equation for the p-spin 
model follow from the h — > limit of Eqs. ( |3.30| ) and ( ^.31[ ). 

APPENDIX D: TREATMENT OF INTEGRALS 



Integrals over time windows spanning the interval [0, t] appear in the dynamic equations 
( |3.30| ) and ( |3.31| ). In the large times limit we approximate them in the way we here describe. 



1. First type of integral 

Integrals of the form 

h(t)= f dt"A(t,f)B(t,t") (D.l) 
J o 

appear, for example, in z(t). The idea is to separate the integration time-interval as 

[0,t] = [0,6] u [6,r] u [r,t] . (D.2) 

If 5 is chosen to be a finite time, all functions can be approximated by A(t, 0) that vanishes 
when t —>■ oo. Since the integration interval is finite, this term can be neglected. In the 
second interval the functions vary in the aging regime and in the third interval they vary in 
the stationary regime. Thus 

f dt"A{t,f)B{t,f) + f dt"A AG (t,t")B AG (t,t") 

JO J5 

+ f dt" \( A ST (t - t") + lim lim A{t, t")) ( B ST {t - t") + lim lim B(t, t")) 

Jt~ \_\ t—t"^oot"^oo J \ t—t"^oot"^oo J 



(D.3) 



50 



We assume that this separation is sharp and that we can neglect the corrections associated 
to mixing of the three regimes. In the third term we replaced A and B in terms of A ST , B ST . 

Since in all cases either A or B is proportional to the response, we can neglect the first 
term using the property ( 4.11|) of the weak long-term memory scenario. We can then replace 



the lower limit of the second integral by and the upper limit of the second integral by t. 
In addition, assuming that B is proportional to the response, 

lim lim B(t,t") = . (D.4) 

t—t"^oo t"^oo 

This yields 

h(t)~ dt"A AG (t,t")B AG (t,t") + ( lim lim A(t,t") / dr' B ST (r') 

JO \t-t"—>oo t"^oo J Jo 

t—t~ -^oo 

dr' A st (t')B st (t') . (D.5) 



o 



2. Second type of integral 

Another type of integrals is: is 

I 2 (t, t') = /' dt"A(t, t")B(t\ t') . (D.6) 
Jt' 

In particular, if B — 1, A(t,t") = R(t,t") and t' = 0, this integrals yields the susceptibility. 
Let us assume that t and t' are far apart; we start as for the first type of integrals, by 
dividing the time interval in three subintervals 

[t',t} = [t',t' + ] U {t' + ,t-} U [t~,t] (D.7) 

and by approximating the functions by their functional form inside each of the intervals: 

I 2 (M')~ f dfA AG (t,t")B(t"-t')+ f dt"A AG {t,f)B AG {t",t') 
Jt' Jt'+ 

+ ^ dt"A(t - t")B AG (t", t') . (D.8) 

In the first term B{t - t") can be replaced by B(t - t") = B ST (t - t") + 
limt_t»^oo lim^^oo B(t, t"). The same applies to A in the last term. All functions vary 
fast in the stationary regime but very slowly in the aging regime. The next assumption 
is that functions in the aging regime that are convoluted with functions in the stationary 
regime, can be considered to be constant and taken out of the integral. That is to say 

h{t, ~ A AG (t, t') f dt" B(t" - t') + f dt"A AG (t, t")B AG (t", t') 
Jt' Jf 

+B AG (t,t')£dt"A(t-t") 

POO ft 

~A AG (t,t') / dr'B(r')+ / dt"A AG (t,t")B AG (t",t') 
Jo Jt' 

+B AG (t,t') / dr'A(r') , (D.9) 
Jo 
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where we used t — t~ — > oo and t' — t' — > oo. Typically, A(r') and -8(1"') are proportional 
to the response function. Using FDT the integrals in the first and third term can then be 
performed in the classical limit or they can be expressed as functions of the correlation in the 
quantum case. The second term instead depends exclusively on the aging dynamic sector. 
All other integrals can be evaluated, in the large-time limit, in a similar way. 
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